Self-interacting Scalar Field Trapped in a DGP Brane: The Dynamical Systems 

Perspective 
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We apply the dynamical systems tools to study the linear dynamics of a self-interacting scalar 
field trapped on a DGP brane. The simplest kinds of self-interaction potentials are investigated: a) 
constant potential, and b) exponential potential. It is shown that the dynamics of DGP models can 
be very rich and complex. One of the most interesting results of this study shows that dynamical 
screening of the scalar field self-interaction potential, occuring within the Minkowski cosmological 
phase of the DGP model and mimetizing 4D phantom behaviour, is an attractor solution for a 
constant self-interaction potential but not for the exponential one. In the latter case gravitational 
screening is not even a critical point of the corresponding autonomous system of ordinary differential 
equations. 
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I. INTRODUCTION 

Since the discovery that our universe can be cur- 
rently undergoing a stage of accelerated expansion 
many phcnomcnological models based either on Einstein 
General Relativity (EGR), or using alternatives like the 
higher dimensional brane world theories Q, have been 
invoked (for a recent review on the subject see reference 
0). The latter ones, being phenomenological in nature, 
are inspired by string theory. 

One of the brane models that have received most at- 
tention in recent years is the so called Dvali-Gabadadze- 
Porrati (DGP) brane world This model describes a 
brane with 4D world- volume, that is embedded into a flat 
5D bulk, and allows for infrared (IR)/large scale modi- 
fications of gravitational laws. A distinctive ingredient 
of the model is the induced Einstein-Hilbert action on 
the brane, that is responsible for the recovery of 4D Ein- 
stein gravity at moderate scales, even if the mechanism 
of this recovery is rather non-trivial Q. The accelera- 
tion of the expansion at late times is explained here as 
a consequence of the leakage of gravity into the bulk at 
large (cosmological) scales, so it is just a 5D geometrical 
effect, unrelated to any kind of misterious "dark energy". 

As with many IR modifications of gravity, there are 
ghosts modes in the spectrum of the theory [7], Hj]. 2 Nev- 



ertheless, studying the dynamics of DGP models contin- 
ues being a very atractive subject of research, ft is due, 
in part, to the very simple geometrical explanation to the 
"dark energy problem", and, in part, to the fact that it 
is one of a very few possible consistent IR modifications 
of gravity that might be ever found. 

The aim of this letter is, precisely, to study the dy- 
namics of a self-interacting scalar field trapped on a DGP 
brane 3 by invoking the dynamical systems tools, which 
have been proved useful to retrieve significant informa- 
tion about the evolution of a huge class of cosmologi- 
cal models (see for instance the book pJp|), The sim- 
plest self-interaction potentials: a) constant potential, 
and b) exponential potential, are investigated. The con- 
stant potential is the simplest one, and, for frozen scalar 
field cf> = cf>o, for the Minkowski cosmological phase, the 
corresponding model coincides with the one of reference 
[9( . The exponential potential represents a common func- 
tional form for self-interaction potentials that can be 
found in higher-order [llj or higher-dimensional theories 
12 1 . These can also arise due to non-perturbative effects 
13j . In addition to the scalar field we consider, also, a 



background fluid trapped on the DGP brane. 

A dynamical study of DGP models with a scalar field 
trapped on the DGP brane has been already undertaken 
in reference [l4[ for an exponential potential, to show 
that crossing of the phantom barrier is indeed possible 
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DGP model; the so called "self-accelerating" branch, or self- 
accelerating cosmological phase 0. 
3 Regarding the so called "normal" branch, or Minkowski cosmo- 
logical phase [ij of DGP model, this is just a particular member 
in the wider class of models of reference [T3| . 
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in DGP cosmology with a single scalar field (see also 
(TH in this regard). However, the authors of that pa- 
per do not study in detail the phase space of the model 
and, in correspondence, they are not able to find critical 
points. Their claim that scaling solutions do not exist 
in a universe with dust on a DGP brane, seems to be 
in contradiction with our results. Also, in reference flr| 
a dynamical systems study of DGP scenarios is under- 
taken for several self-interaction potentials. However the 
author of the latter reference considers only a scalar field 
trapped on the brane (no other kind of matter is trap ped 
on it). Besides, phase space variables used in [16[ are 
different from the ones used in the present study. 

In a sense, the present investigation is a continuation 
of the one reported in paper [17] . to include higher- 
dimensional behaviour (in the present case dictated by 
the DGP dynamics). In consequence, for the exponen- 
tial self-interaction potential, our results will include the 
ones reported in [171 ] as a particular case. 

Through the paper we use natural units (87rG = 
87r/mp ; = h = c = 1). 



following choice of dynamical variables: 4 

1 w 



X v/6Q' V ~ V3 Q ' *T' [ ) 

and, also, we introduce the time variable r = J Hdt. 
The following autonomous system of ordinary differential 
equations is obtained: 



x = -3x-^ (d^lnV) y 2 z + 

^x( 1 (l-x 2 ~y 2 ) + 2x 2 ) 1 
y ' = Vl ^ lnl/ ) xyz+^y( 1 (l-x 2 -y 2 ) + 2x 2 ), 
Z ' = \ Z (^)(7(l-* 2 -y 2 ) + 2z 2 ), (4) 

where a prime denotes derivative with respect to t. In 
terms of the dynamical variable z, equation @ can be 
rewritten as: 



II. THE MODEL 



z z = l± 



1 



(5) 



As already mentioned, we will be concerned here with 
a DGP brane model where self-interacting scalar field 
matter is trapped on the DGP brane. The field equations 
are the following: 



i 



(p+ p + V {fl>)), 



p=-3 7 Hp, <t>=-3H<t>-d (f ,V, 



(1) 



where p is the energy density of the background 
barotropic fluid (7 is the barotropic index) , 4> is the scalar 
field trapped in the DGP brane, V its self-interaction po- 
tential, and we have used the following definition: 



Q% = H 2 ± — H. 



(2) 



with r c being the so called crossover scale. Depending on 
the choice of the signs " +" or " -" in @ , there are two pos- 
sible branches of the DGP model, corresponding to two 
possible embeddings of the DGP brane in the Minkowski 
bulk. The choice " +" is for the so called Minkowski cos- 
mological phase of DGP that is free of ghosts, while the 
"-" choice is for the self-accelerating cosmological phase. 
The Minkowski phase of the present model belongs in the 
wider class of models of reference 15]. For <j) — const, the 
self-accelerating phase coincides with the model of 



For the Minkowski phase, since < H < 00 (we consider 
just non-contracting universes), then 1 < z < 00. The 
case —00 < z < — 1 corresponds to the time reversal 
of the latter situation. For the self-accelerating phase, 
— 00 < z 2 < 1, but since we want real valued z only, 
then < z 2 < l. 5 As before, the case — 1 < z < 
represents time reversal of the case < z < 1 that will 
be investigated here. Both branches share the common 
subset (x, y, z — 1), which corresponds to the formal limit 
r c — > 00 (see equation ([2])), i.e., this represents just the 
standard 4D behaviour of (4D) Einstein-Hilbcrt theory 
oupled to a self-interacting scalar field. 

In terms of the above introduced dynamical variables 
x,y,z, the "effective" dimensionless density parameters 
tt = p/3Q 2 , and fi = p^/ZQ 2 {p<p = <j) 2 /2 + V) can be 
expressed in the following way: 



i-x 2 - y 2 , n 4 



*-* 2 + y\ (6) 



where, as customary, Q = p/3H 2 , Vl^ 



h /3H 2 . Since 



both tt > 0, > 0, and z z > 0, then < x z 



y 



< 1. 



This means the phase space for the dynamical system 
driven by the evolution equations ([3]) , for the Minkowski 
cosmological phase (the ghost- free " +" branch) , is given 
by the unbounded region: 



III. THE PHASE SPACE 

Our aim is to write the system ([T]) as an autonomous 
system of ordinary differential equations. We make the 



4 In what follows, for brevity, we avoid writing the "±" sign. 

5 In fact, fitting SN observations requires H > in order to 
achieve late time acceleration (see, for instance, reference Q and 
references therein). This means that z has to be real- valued. 
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A. The Constant Potential V = Vo 



^+ = {{x,y,z): 0<x 2 + y 2 < 1, ze [l,oo[}, (7) 

while, for the self-accelerating cosmological phase (the 
pathological "-" branch), it is given by the complemen- 
tary (non-compact) region: 



= {(x,y,z): 0<x 2 +y 2 <l, ze]0,l}}- (8) 



The autonomous system of ordinary differential equa- 
tions (HJ) simplifies to: 



x' = -3x + -x(7(l - x 2 - y 2 ) + 2x 2 ), 

V' = §2/(7(1 - x 2 - y 2 ) + 2a; 2 ), 

z' = \z(i r ±){ 1 {l-x 2 -y 2 ) + 2x 2 ). (11) 



Notice that the points belonging in the set (a;, y, 0) can 
not be included since, in this case {z = Q = 0) 
the variables x and y are undefined. In consequence, the 
self-accelerating solution H — l/r c has to be estudied 
separately, i. e., with a different choice of phase space 
variables (see section |V)) . 

Take a look at (j4|). These equations hide the differ- 
ences between the "+" and "-" branches of the DGP. 
Nevertheless, the differences indeed exist. These are 
encoded in the different definition of the corresponding 
phase spaces (see the definitions for ty+ and above) 
which, from now on, we will call as Minkowski and 
self-accelerating (4 r _) phases respectively. Recall that 
points in the phase plane {{x,y, 1) : < x 2 + y 2 < 1}, 
where both and intersect, are related to standard 
4D behaviour (formal limit r c — > oo). 

To finalize this section we write another useful mag- 
nitudes of observational interest in terms of the dynam- 
ical variables ||3J). These are, the deceleration parameter 
q= -(l + H/H 2 ): 



3z 2 



z 2 + l 



( 1 (l~x 2 -y 2 ) + 2x 2 ), 



(9) 



and the equation of state (EOS) parameter w$ 
2V)/(<P + 2V): 



2 2 

xr — y z 
x 2 + y 2 



(10) 



IV. THE CRITICAL POINTS 

We study here two concrete yet generic cases. First we 
consider a constant potential V — Vq, and then we study 
the exponential potential V — Vq exp (— a<f>). For more 
general classes of potentials we have to rely on a quite 
different approach (see, for instance reference [HI). This 
will be the subject of future research. 



Before finding the critical points of this system of equa- 
tions and investigating their stability, a straightforward 
inspection of the differential equations (fTTj) reveals that 
the first two equations (differential equations for x and y 
variables respectively) are decoupled from the third one. 
This means that there can be found critical points with 
x = 0, y = ±1 and arbitrary z-s to form critical sub- 
manifolds of particular interest (see below on the critical 
submanifolds M±). These critical points are associated 
with the inflationary (q — —1) DGP-Fricdmann equation 
(notice that, in this case, Q — 0, = z 2 ): 



(12) 



H 2 ±i H =m. 

r c 3 



For the Minkowski cosmological phase ("+" sign in ljl2l ). 
this solution can be associated with 4D phantom-like be- 
haviour, that is generated through dinamical (gravita- 
tional) screening of the potential energy V(4>). 6 As we 
will see in the next sections, the latter solution is an at- 
tractor only for the constant potential V — Vq. For the 
exponential potential it is not even a critical point of the 
corresponding autonomous system of ordinary differen- 
tial equations. 

The critical points of the autonomous system of equa- 
tion (fTTj) are summarized in table HI while table HT1 shows 
the eigenvalues of the corresponding Jacobian matrices. 
Notice, from (|11[) . that the dynamical equations are in- 
variant under the change of sign z — > — z, in consequence 
we have not included the points with z = — 1 in our anal- 
ysis. 

There are 5 critical points Pi, P 2 ± , and two criti- 
cal submanifolds M± = (0,±l,z) (z e]0, 1] for the self- 
accelerating phase, while, z £ [l,oo[ for the Minkowski 
phase). It is obvious that critical points P 3 ± belong in 
M±. 

Notice that for critical points P\, P, 



Pg , SHUT 



2 ' 



z = 1 (formal limit r c — > oo), then Q = Q, = Q,p, 
i.e, these points correspond to standard 4D behaviour. 
Otherwise, critical points attached to the phase plane 
{(x,y, 1) : < x 2 + y 2 < 1} are the well known critical 



The present case (V = Vo) coincides with the case studied in [jj. 
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TABLE I: Properties of the critical points for the autonomous system 



p 


X 


y 


z 


Existence 


Op 


U)<f> 


q 


Pi 








1 


Always (V7 £ [0,2]) 





undefined 


-1 + 321 


p ± 


±1 





1 


•>i 


1 


1 


2 


p ± 





±1 


1 


w 


1 


-1 


-1 


M± 





±1 


z e]o, oo [ 


)) 


1 


-1 


-1 



TABLE II: Eigenvalues for the critical points in table U 



Pi 


X 


y 


z 


Ai 


A 2 


A3 


Pi 








1 


-3(2- 7 )/2 


3 7 /2 


3 7 /2 


p 2 ± 


±1 





1 


3(2-7) 


3 


3 


p? 





±1 


1 


-3 


-3 7 





M± 





±1 


z 


-3 


-3 7 






points for, basically, standard 4D Einstein-Hilbert grav- 
ity with a minimally coupled self-interacting scalar field 
with constant self-interaction potential. These points ex- 
ist for both branches of the DGP model. 

The matter-dominated solution Pi (f2 = fi = 1) is 
a saddle point in $±. It corresponds to an accelerat- 
ing solution for 7 < 2/3 (decelerating otherwise, i.e, for 
2/3 < 7 < 2). Stiff-matter dominated solutions P 2 ± 
(Clcj, = CIq = 1, cj0 = 1) represent source (repellor) crit- 
ical points in phase space, and are always decelerating 

(<z = 2). 

The scalar field dominated inflationary phases (f2^ = 
il0 = 1, W0 = — 1, q = — 1) corresponding to points in 
the submanifolds M± in table U (as already said, these 
include the critical points P^), represent non- hyperbolic 
critical points, since one of the eigenvalues of the corre- 
sponding Jacobian matrixes vanishes (see table [TTJ) . In 
this case the only thing we can state with certainty, on 
the basis of straightforward analisys of the autonomous 
system of equations (fTTj) . is that trajectories in phase 
space, originating in one of the repellor points P 2 ± , de- 
pending on the phase considered: the Minkowski phase 
or the self-accelerating one, and on the initial condi- 
tions, will inevitably approach one or several of the points 
(0, ±l,Zjo) € M±. Otherwise, points in the segments 
(0, ±l,z), that are parallel to the z-axis, can be seen 
as attractors by the corresponding phase space trajec- 
tories. Actually, the stable subspace of the dynamical 
system near of the critical points in the submani- 

folds (0, ±1, z), is spanned by the eigenvectors: 





(A 
















Co) 







i.e., the stable subspace is the plane (%,y) intersecting 
the z-axis at a given z = z^. The central subspaces co- 
incide with the open segments (0, ±1,2 G]0, 1]) for the 
self-accelerating phase, while, for the Minkowski phase 
these coincide with the semi-infinite segments (0, ±1, z £ 
[1, oo[). These segments are parallel to the z-axis and 
lie on the boundary of 1 $>± (9^±). Since, in this case, 
there is no unstable subspace related to M±, and the 
centre manifold belongs in the boundary of the phase 
space, this means that trajectories in phase space that 
evolve from the stiff-matter repellor £ (x,y,l), will 
inevitably approach one or several points (depending on 
the cosmological phase being considered) in M± , into the 
future (r — > 00). As we will see quite soon, it is a fact 
that, only for the Minkowski phase the points in 
M± are attractor critical points, so that, the correspond- 
ing submanifolds are attractor segments. These critical 
subspaces are an important result, since, as we will see 
in the next subsection, they exist only for the constant 
potential case, meaning that the associated cosmological 
solutions are generic only in this particular case. 

Additional and, perhaps, more precise information can 
be extracted from the phase portraits. Actually, in the 
figure HJ phase trajectories are shown for the constant 
potential V = Vo, for different sets of initial conditions. 
The upper figure is for the self-accelerating phase 
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while the figure in the middle is for the Minkowski phase 
The lower figure shows the flow of the autonomous 
system (fTTj) . From these figures several features are ap- 
parent: 

• Phase trajectories in , F_ (upper figure in fig. [1]) 
originate from the source critical points P 2 ± (°nly 
the point is shown), correponding to the stan- 
dard 4D kinetic energy dominated (stiff-matter) so- 
lution, and (asymptotically) approach to the point 
(0,1,0) that has been removed from the phase 
space since phase space variables x and y blow up 
at the phase plane (x,y,0). The dynamics in the 
neigbourhood of this point has to be investigated in 
terms of different phase space variables. Therefore, 
for the self-accelerating cosmological phase "J- the 
points (0, ±1, Zi) £ M± are not even critical points 
of the autonomous system (TTT]). 

• Phase trajectories in ^ + (center figure in fig. [T]) 
originate from the 4D stiff-matter solution (unsta- 
ble node P£ in tab. [TJ) and end up at the infla- 
tionary points (0, 1, ZQi) £ M+, where the different 
z 0i -s are associated with the different initial condi- 
tions. Otherwise, points in M + are seen as attrac- 
tor points by the different phase space " observers" , 
moving along different phase trajectories that orig- 
inate at P 2 + . These scalar field dominated critical 
points correspond to inflationary solutions of the 
DGP-Friedmann equation: 



That the phase trajectories in \t + (middle figure in 
fig. [T|), leave the phase plane (x, y, 1) and probe the bulk 
of the phase space >F+, is a nice illustration of the fact 
that 4D phantom-like behaviour arising from dynamical 
screening of the brane cosmological constant Vb [9( , is a 
phenomenom of 5D nature. The new feature revealed 
by the present dynamical systems analysis relies on the 
fact that this type of behaviour - gravitational screen- 
ing - is quite independent of the initial conditions, since 
the corresponding semi-infinite segments M± £ are 
stable submanifolds. 




FIG. 1: Trajectories in phase space for different sets of initial 
conditions for the constant potential V = Vo- The upper 
figure and the figure at the center are for the self-accelerating 
and Minkowski >]/+ phases of the DGP model respectively. 
The figure at the botton shows the flow in time r. Notice 
that trajectories in asymptotically approach to the self- 
accelerating solution H — l/r c . Meanwhile, the trajectories 
in ^+ approach points in the segment (0, 1, z). 



B. The exponential Potential V = Voexp(— aq 



I [3 2 ^ 2 2 2 

The autonomous system of ordinary differencial equa- x ~ ^ x ^~ \J V z + ~^ x {l{^ ~ x ~ V ) + 2x ), 

tions (TJ| takes the following form: 



7 This argument has been suggested to us by Roy Maartens. 2 



3 3 

- a xyz + -y(7(l - x 2 - y 2 ) + 2x 2 ), 
z' = ~z(^)(j(l-x 2 - y 2 ) + 2x 2 ). (13) 
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TABLE III: Critical points for the autonomous system l|13[l. 



Pi 


X 


y 


2 


Existence 






9 




Pi 








1 


All a, 7 





undefined 


-1 + 


3l 
2 


P 2 ,3 


±1 





1 


All a, 7 


1 


1 


2 




Pi 


a 




1 


a 2 < 6 


1 


<£_! 
3 ± 


a 2 
2 


1 


P 5 


[a i 

V2o 




1 


a 2 > 37 




7-1 


-1 + 


2 



Before going any further into the phase space analysis 
of (|13p . let us note that, contrary to the case with the 
autonomous system of diferential equations (jllj) . thanks 
to the term with the factor d t f,V(4>) = —a in the right- 
hand-side of the first couple of equations, the equations 
in (fl~3"|) form a coupled system of differential equations. 
The inmediate consequence is that points in the subsets 
(0, ±l,z) £ d*f?± are not even critical points of the au- 
tonomous system (p~3|) . 

The critical points of (| 13[) are summarized in table HTT1 
We have not included the points with z = — 1, since these 
coincide with those with z = 1 under the change of sign 
of the exponent a — ► —a. Recall that the critical points 
with z = 1 correspond to the formal limit r c — > 00, i. 
e., these represent standard 4D behaviour. A complete 
study of the critical points in table IIIII can be found in 
the well known reference [l7j . The main results can be 
summarized as follows: 

• a 2 < 37; The kinetic-dominated solutions (points 
P 2i 3 in table are unstable nodes (sources). The 
matter-dominated solution (point Pi) is a saddle 
point. The scalar field dominated solution (point 
P4) is the late-time attractor. 

• 37 < a 2 < 6; The kinetic-dominated solutions P2,z 
are unstable nodes. The scalar field dominated so- 
lution P4 is a saddle. The matter-scaling solution 
(point P5 in table IIII| ) is a stable node/spiral. 

• a 2 > 6; The kinetic-dominated solution can be ei- 
ther a source or a saddle. The matter-dominated 
solution is a saddle point. The matter-scaling solu- 
tion P5 is a stable spiral. 

The main difference with the case for the constant po- 
tential V = Vq of subsection IIV Al is in the absense of the 
critical submanifolds M± that can be associated with 4D 
phantom-like behaviour originated from 5D dynamical 
screening of the potential energy of the scalar field. Al- 
though 5D behavior related to dynamical screening does 
actually arise in the present case (as can be seen from the 
figure the phase trajectories leave the plane (x, y, 1) 
and probe the bulk of the 3D phase space as for the 
constant potential), the corresponding DGP-Friedmann 



behaviour dictated by the equation 

H 2 + -H = \ 

is not even a critical point, so that 5D gravitational 
screening mimetizing 4D phantom behaviour is not a 
generic solution in this case. Actually, it is evident from 
the middle figure in Fig. [2l that the trajectories in phase 
space *f? + (the Minkowski cosmological phase) do not 
even touch the segment (0, 1, z) £ d^+. In fact, for the 
chosen initial conditions, these seem to be repelled from 
that segment. 

Another interesting conclusion that can be made after 
the phase space a naly sis in this subsection is that, the 
claim in reference [14| that scaling solutions do not exist 
in a DGP brane filled with dust, in general, is not correct. 
For a 2 > 37, matter scaling solutions indeed exist, they 
are associated with critical points in phase space (point 
P5 in table HIT) ) . and can be, even, attractor solutions. 



V. SELF-ACCELERATING SOLUTION H = l/r c 

As already noted from the former analysis, for the self- 
accelerating phase of the DGP brane model, independent 
of the initial conditions, trajectories in phase space ap- 
proach to the self-accelerating solution with phase space 
coordinate z = =>• H = l/r c . However, the points 
with z — are to be removed from the phase space 
since, otherwise, the coordinates x and y in ((3]) are unde- 
fined. Therefore, a new choice of phase space coordinates 
is mandatory to find the stability properties of the self- 
accelerating solution in phase space. Let us start with the 
DGP-Friedmann equation for the self-accelerating cos- 
mological phase, written in the following general form 



H 2 



H 



(14) 



where p is now the energy density of the total matter 
content trapped on the DGP brane. If new phase space 
variables are defined: 



r,iT 



V3P' 



(15) 
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FIG. 2: Phase trajectories for different sets of initial condi- 
tions for the exponential potential V = Vbexp(— ac/>). As be- 
fore the upper figure and the figure at the center are for the 
self-accelerating ^_ and Minkowski ^+ phases of the DGP 
model respectively. The flow of the autonomous system in 
time r is shown in the figure at the botton. Unlike the case 
for the constant potential in subsection IIV Al phase trajecto- 
ries in *&+ are repelled by the segment (0, 1, z) (middle figure). 



then, the following Friedmann constraint arises: 

l=x + if =*> y = Vl-x, (16) 

so that, to have real valued y, the variable x is con- 
strained to the interval x E [0,1]. This, in turn, yields 
that y E [0, 1] (the case y E [—1,0] is related with a con- 
tracting phase of the cosmic evolution) . In consequence 



the phase space is = {(x,y) : x E [0,1], y E [0,1]}. 
Due to the Friedmann constraint p^|) . we are left with 
only one independent autonomous differential equation: 

z' = 37z(^£). (17) 

The corresponding critical points are the standard (4D) 
Friedmann solution Pi = (0,1) => H 2 = p/3 and the 
self-accelerating (5D) solution P 2 = (1,0) => H = l/r c . 
To study the stability we must perturb (|17j) in the neigh- 
bourhood of the critical points: x — > x\ + e. Consider 
linear perturbations, i. e., we neglet terms ~ C(e 2 ), then 

e ' = 3 ^' + e >(}^)- < 18 > 

so that, for the Friedmann solution (point Pi), 

e ' = ^ e + 0( e 2 ) => e (r)=e(0)e^, (19) 
while, for the self-accelerating phase: 

e' = -3 7 e + 0{( 2 ) => e(r) = e(0) e^ 7 ". (20) 

From the analysis of (|19p and ([20| it is evident that the 
standard Friedmann solution is unstable, while the self- 
accelerating one is stable against small (linear) pertur- 
bations. Otherwise, in the languaje of the dynamical 
systems analysis, the latter solution is an attractor in 
phase space. This solution is, precisely, the one corre- 
sponding to the point (x,y 7 z) — (0,1,0) in the former 
section, that has been removed from the correponding 
phase space. Therefore, the reason why phase trajecto- 
ries in 'F- of the former section (for both constant and 
exponential potentials) approach to that point, is that, 
in adequate phase space coordinates, the self-accelerating 
solution H — l/r c is an attractor critical point. 

VI. RESULTS AND DISCUSSION 

From the analysis in the former sections, the following 
important results can be summarized: 

• For the constant potential the first 5 critical 
points (Pi, P^ and P 3 ± ) coincide with those in 
standard 4D Einstein-Hilbert theory coupled to 
a self-interacting scalar field with constant self- 
interaction potential. However, in principle, the 
phase trajectories are not constrained to the phase 
plane (x, y, 1) as in standard 4D gravity coupled to 
a self-interacting scalar field, these probe the bulk 
of the 3D phase space. For solutions in >F_, the 
(standard 4D) stiff-matter solution (critical points 
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P 2 ± in table H| are sources (unstable nodes), while 
the 4D matter-dominated solution (point Pi) is a 
saddle point. Phase trajectories (asymptotically) 
approach to the point (0, 1, 0) that has been re- 
moved from the phase space to avoid singularities 
related to the choice of phase space variables x, y 
and z. This point, when studied separatelly - with 
a different choice of phase space coordinates - cor- 
responds to the self-accelerating solution H = l/r c , 
which is stable against small (linear) perturbations, 
i. e., it is an attractor solution for phase trajecto- 
ries in 

• There are two critical segments M± = (0, ±l,z) 
that belong in the boundary of the Minkowski cos- 
mological phase 9^+ and exist only for the con- 
stant potential case. Critical points in these sub- 
sets of the phase space correspond to inflationary 
solutions and are non- hyperbolic. Only stable and 
centre submanifolds can be attached to them. A 
strightforward analysis of the autonomous system 
of ordinary differential equations (fTT|) . shows that 
the stability of these critical subspaces depends 
on the way phase space trajectories approach to 
them. Phase portraits, however, reveal that points 
in these segments are seen by the corresponding 
phase trajectories as sinks or attractor points. 

• In general, the dynamical behaviour of the DGP 
model with exponential self-interaction potential, 
does not differ from the standard behaviour of 
4D Einstein-Hilbert gravity coupled to a self- 
interacting scalar field (see reference Qj}), but for 
the fact that trajectories in phase space, inevitably, 
leave the phase plane (x, y, 1) and probe the bulk 
of the 3D phase space due to the 5D nature of the 
DGP model. In this case, however, there are not 
critical points that can be associated with (5D) 
gravitational screening of the scalar field energy 
density. 

Perhaps, the most interesting finding of the present 
investigation can be associated with the existence of the 
critical submanifolds M± — (0, ±l,z £ [l,oo[), that are 
distinctive of the constant potential case only. The crit- 
ical points in M± are associated with dynamical screen- 
ing of the cosmological constant Vq and, since these are 
attractors in the fate of the associated cosmologi- 
cal evolution is quite generic and independent of the ini- 
tial conditions. Therefore, the existence of M± serves as 
a nice illustration of dynamical screening of the cosmo- 
logical constant (a phenomenon! of 5D nature) in phase 
space. The fact that, for the exponential potential case, 
the semi-infinite segments (0,±,z) £ are not even 

critical subspaces, came as a surprise. The inmediate 
consequence is that, for the exponential potential, the 
gravitational screening is not as generic as for the con- 
stant potential. 

The latter conclusion is nicely illustrated in the figure 
|31 where phase trajectories have been projected onto a 




FIG. 3: Projection of phase trajectories onto a given plane 
(x,y,zo) for the Minkowski phase of the DGP model, for 
the constant potential (upper figure), and for the exponen- 
tial potential (lower figure), respectively. We chose an arbi- 
trary zq = 2. Notice that, while the point (x,y) = (0,1) 
is an attractor for V — Vo (upper figure), for the case 
V — Vbexp(— acf>) (lower figure), it is not. 



given phase plane (x,y,Zo)- It is seen that, for the case 
V{4>) = Vq, trajectories in phase space converge towards 
the point (0,1) (recall that we fixed z = zq), while for 
V{4>) — Vq exp(—a(p) these trajectories are repelled by 
this point. 

The fact that phase trajectories probe the bulk of the 
3D phase space means that, depending on the initial 
conditions, the cosmological dynamics of DGP models 
can be very rich and complicated. A clear illustration of 
this assertion is associated, precisely, with the occurrence 
of the attractor (self-accelerating) solution H = 1 /r c in 
and of the critical (also attractor) subspaces M± £ 
d^+, for the constant potential case. 



VII. CONCLUSION 

The dynamics of DGP models can be very rich and 
interesting. A dynamical systems approach to the sub- 
ject reveals that trajectories in phase space scape from 
the plane (x, y, 1), that is associated with standard 
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4D Einstein-Hilbert theory coupled to a self-interacting 
scalar field. This feature is associated with the 5D na- 
ture of the DGP model and originates a phenomenom 
called as dynamical screening, that mimetizes 4D phan- 
tom behaviour. While the latter phenomenom is quite 
generic and independent of the initial conditions if the 
self-interaction potential is a constant (V = Vq), when 
V — Vq exp(— a<f>), the gravitational screening indeed oc- 
curs, but it is not a generic phenomenom. 

Another interesting result can be associated with the 
existence of matter-scaling solutions that can have, even, 
attractor nature. This result is in contradiction with 
the claim in [l4[ that scaling solutions do not exist in 
DGP universes with dust. In the case of the exponential 
potential, the existence of scaling solutions is expected 
since, the corresponding DGP model is a generalization of 
the Einstein-Hilbert theory coupled to a (self-interacting) 
scalar field with V — Vo cxp(— acf) (the dynamics of this 
model has been investigated in reference [17|). to include 
5D effects originated from induced gravity on the brane. 
Then, since in the formal 4D limit r c — > oo the results 
of reference [llj are to be recovered, the critical points 
found in that reference - including the scaling solution 



- have to be, also, critical points of the corresponding 
DGP model. 

It can be of insterest to investigate the present DGP 
scenario for arbitrary self-interaction potentials, to show 
that independence of the gravitational screening of the 
initial conditions is distinctive only of the constant po- 
tential case. This task would entail a different approach 
than the one undertaken in this paper, so that we leave 
it for future work. 
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